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Abstract 

It has been suggested that the high symmetries in the Schrodinger equation with the Coulomb or 
harmonic oscihator potentials may remain in the corresponding relativistic Dirac equation. If the 
principle is correct, in the Dirac equation the potential should have a form as V{t) where V{r) 
is for hydrogen atom and Kr^ for harmonic oscillator. However, in the case of hydrogen atom, 
by this combination the spin-orbit coupling term would not exist and it is inconsistent with the 
observational spectra of hydrogen atom, so that the symmetry of S0(4) must reduce into SU(2). 
The governing mechanisms QED and QCD which induce potential are vector-like theories, so at 
the leading order only vector potential exists. However, the higher order effects may cause a scalar 
fraction. In this work, we show that for QED, the symmetry restoration is very small and some 
discussions on the symmetry breaking are made. At the end, we briefly discuss the QCD case and 
indicate that the situation for QCD is much more complicated and interesting. 
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I. INTRODUCTION 



The Schrodinger equation is the basis of non-relativistic quantum mechanics. It success- 
fully explains the energy spectra of hydrogen atom and harmonic oscillator, moreover as the 
orbit-spin coupling is taken into account, the fine structure of hydrogen atom is also calcu- 
lated. All the predictions are consistent with measurements of high precision. In a complete 
theoretical framework which respects the Lorentz covariance, the relativistic Klein-Gorden 
equation for spin-zero particles and Dirac equation for spin 1/2 particles substitute the 
Schrodinger equation to describe physics. Surely, under the non-relativistic approximation, 
they can reduce back to the Schrodinger equation and meanwhile the orbit-spin coupling 
and relativistic correction terms which were added into the non-relativistic Hamiltonian by 
hand, automatically appear. Nowadays, as the detection technique and instrumentation are 
ceaselessly improved, small deviation between theoretical predictions and observational data 
may be exposed, such as the Lamb shiftQ, fl , and it enables one to investigate the tiny 
difference caused by various theoretical models. 

No wonder, as is well known, symmetry is extremely important in all fields of physics. 
Generally, the Schrodinger equation for hydrogen atom which is of a Coulomb-type potential 
of N dimensions possesses an S0(N+1) symmetry. It is not a geometric symmetry, but a 
dynamical one. The SO (4) group in three dimensions can be decomposed into a direct 
product of two independent SU(2) groups which correspond to spin and orbital angular 
momenta respectively 4|. Recently, Zhang, Fu and Chenj^ indicate that the Schrodinger 
equation with a central force fields has an explicit SO (4) symmetry, but generally the Dirac 
equation with the same potential does not possess the SO (4) symmetry. The harmonic 
oscillator has the similar situation. The Schrodinger equation has an SU(3) symmetry, but 
the corresponding Dirac equation does not. Ginocchiop, 3] investigated this problem and 
obtained a modified Dirac equation which possesses the SU(3) symmetry and can also reduce 
into the original Schrodinger equation under non-relativistic approximation. 

2 

In the central force field the non-relativistic Hamiltonian is H = + V{r), turning into 
the operator form, one only needs to replace p into —iKV and does not change the potential 
V{r) in the configuration representation, the Schrodinger equation of the central force would 
retain the symmetry of V{r) for hydrogen and nr'^ for harmonic oscillator). By contrast. 
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in the relativistic theory, the Hamiltonian of a free particle should he H = y/p^ + m^, thus 
if one writes it into an operator form, the square root a/p^ + needs to be replaced by 
a ■ p + m(3. However, if there was a potential V{r) in the non- relativistic Hamiltonian, then 
how to write a corresponding term in the Hamiltonian when the Dirac equation replaces the 
Schrodinger equation, would be a question. In other words, one can ask whether he needs to 

n 

change its form. The authors of Ref . [5| showed that if we do not make a proper changes on 
the form of V{r), the corresponding Dirac equation with the central force filed fails to possess 
the original S0(4) symmetry at all. Alternatively, if we re-write it as V{r){l + /3)/2, the 
SO (4) symmetry would be maintained. Generally, The non-relativistic theory is a natural 
approximation of the relativistic theory, thus one may suppose that the symmetry degree in 
the relativistic theory should be higher than the non-relativistic one. Therefore, there is a 
compelling reason to believe that the Dirac equation of hydrogen atom should respect the 
SO (4) symmetry. 

However, the physics does not promise the symmetry argument. As the authors of Ref.j^ 
indicate, when the Hamiltonian in the Dirac equation is of the form : a-p+(3M+^{l+(3)V , as 
it reduces into the Schrodinger equation, the spin-orbit coupling term does not exist, namely 
the 2Pi/2 and 2P3/2 states should be degenerate, i.e. the fine structure of hydrogen would 
not be experimentally observed. That conclusion is definitely contrary to the experimental 
measurement. 

Let us investigate the origin of the scalar and vector potentials in the Dirac equations 
based on the quantum field theory (QFT). Considering all possible Lorentz structures of the 
interactions which conserve parity, there are five types of interactions between two charged 
fermions Q]: 1 (g) 1 ~ Vs{r); 75 ® 75 ~ 0; 7^ ® 7^ ~ VV(r); 7''75 ® 7^75 ~ Si ■ S2VA(r) 
and cr^jy ® a'^'^ ~ Si ■ S2VT{r). As is believed, the underlying theory which is responsible for 
electromagnetic processes is QED which is a vector-like gauge theory. At the leading order 
only 7^ ®7/i exists, which induces merely the vector potential Vv{r) by exchanging a virtual 
photon. However, the higher order effects, i.e. multi-photon exchanges can cause scalar and 
tensor interactions. If we ignore the part VA{r) and Vrir) which are related to the hyper- fine 
structure in hydrogen atom case and tiny, only the scalar potential should emerge. 

The observation raises a serious problem that the high symmetry in the Schrodinger 
equation which is a non-relativistic approximation of the Dirac equation might be superficial. 
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When the spin-orbit couphng is taken into account, the symmetry would disappear. In other 

2 

words, in the Schrodinger equation with only the Coulomb's potential the apparent 
SO (4) symmetry does not survive in the relativistic Dirac equation. Interesting issue is how 
much the symmetry can be restored. In fact, in some literatures which is aiming to study 
phenomenology, one writes the potential as {b + aj3)V{r) with a + b = 1, and lets a be a free 
parameter to be fixed by experimental data. 
With the form provided in Ref. 
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and by our above discussion, we write it as if = a ■ p + [3M + {b + a[3)V , then following 
Greiner's calculation which includes all possible combinations of a and 6, we can use the 
hydrogen atom spectrum which has been measured with high precision to determine a and b 
and see to what degree the S0(4) symmetry can be restored or never. We find that for the 
hydrogen atom, the fraction of the scalar potential is only at order of 0.1%, namely in the 
QED case (responsible for the hydrogen atom) the symmetry breaking is almost complete, 
i.e. ^0(4) — ^ SU{2). 

Based on the quantum field theory, the global and local symmetries which reside in 
the Lagrangian are real, instead, the symmetries in the Schrodinger equation are indeed 
superficial. On other aspect, the higher effects can partly restore the symmetry in the 
potential, in fact, the information can help us understanding the effects of higher orders in 
QED and even in QCD where higher orders and non-perturbative effects may play a crucial 
role. 

After this long introduction, we will present our numerical results on the hydrogen atom 
spectrum and determine the fraction of scalar potential in section II. Then we analyze the 
origin of the scalar potential based on the QED theory and discuss the breakdown of SO (4). 
The last section is devoted to our conclusion and meanwhile we briefly discuss the QCD 
case and indicate that for QCD the situation is much more complicated and the symmetry 
analysis might help us to get a better understanding of the non-perturbative QCD effects 
and higher order contributions. 
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II. DIRAC EQUATION FOR HYDROGEN ATOM 



In order to further explore the symmetry problem, we investigate the symmetry issue of 
the most familiar hydrogen atom which has been thoroughly studied from both theoretical 
and experimental aspects. In the case of hydrogen atom, the interaction between the electron 
and proton is purely electromagnetic, so that one can determine to what degree the SO (4) 
symmetry is broken. 

The Dirac equation with an unbroken SO (4) symmetry has the form 



[ca-p + mc^/3+^-^V{r)- Elip = 0, (1) 

and for hydrogen atom V{r) = — 

In this expression the vector and scalar potentials coexist with the same weight. Its 
solution is explicitly given in Ref. ^, 

[ca-p + (3{m^ + V2) -{E- V^)]^ = 0, (2) 

where Vi and V2 are vector and scalar potentials respectively. If one chooses Vi = V2 = 
—hc-^ = — 1^, he can obtain EqO The spectrum is 

E = m[l--^^A (3) 



a 



4n2 



where the natural unit system c = h = 1 is adopted and the principal quantum number is 
n = n' + j + |. 

From Eq.(l3]) we can see that the energy spectrum is determined only by the principal 
quantum n i.e. the states with various j— values are degenerate and there would be no fine- 
structures, and it is the picture described by the Schrodinger equation whose potential does 
not involve the spin-orbit couplings. 

Generally the Dirac equation for hydrogen atom as only the vector potential being in- 
volved, takes the form which is given in ordinary textbooks of quantum mechanics: 

e2 

[ca-p + mc^p- — -E]^p = 0. (4) 
The corresponding eigen-energies are 

^ = '"[1 + tj!^!-'. (5) 



where s = ^J {j -\- \Y — and the energies are exphcitly related to quantum number j. 

We present the eigen-energies corresponding to various j and / when n = 1,2 in the 
following table, where Ei and E2 are calculated in terms of Eqs.Q and ([3]) respectively and 

TABLE I: 



n 


n' 


j 


K 


/ 


structure 


£;i(eV) 


^2(eV) 


1 





1 
2 


1 





Isi 


13.605433412 


13.605071149 


2 





3 
2 


2 


1 


2p3 


3.400994390 


3.400971749 


2 


1 


1 

2 


-1 


1 


2pi 


3.401039674 


3.400971749 


2 


1 


1 
2 


1 





2si 


3.401039674 


3.400971749 



the mass of electron and fine-structure constant are set according to Ref. [15|. Since the 
spectrum of hydrogen has been experimentally measured with very high accuracy, we keep 
more significant figures in our calculation results accordingly. It is shown later that some 
small distinctions between theory and data only appear at the last few digits of the number. 



III. SYMMETRY IS RESTORED PARTLY OR NEVER 

As discussed above, let us rewrite the Dirac equation as 

[ca-p + mc^(3 + hV{r) + a/5V(r) - E]^ = 0, (6) 

with a constraint a + b = 1 which ensures this equation to reduce back to the familiar 
Schrodinger equation, but now we suppose a and b to be free parameters (in fact, only a or 
b is free by the constraint). 

It is easy to obtain the eigen-energy of hydrogen according to Ref. ^ as 

— (1 — a)aa^ (1 — a)aa'^ 2 

= ""Vl - a)2a2 + fn - 7 - K + s)2 + n _ .^2^,2 ^ _ , _ i j_ ,^2- 



aya^ + [n - J - ^ + s)' (1 - a^a^ + [n - j - ^ + s) 
ci^a^ — in — j — -\- sY 1 



2 



}, (7) 



(1 - aYa^ + {n- j -\ + sy 

Generally, one can determine the parameter a by fitting the data of hydrogen atom. It 
is noted that all the measured data of hydrogen atoms involve the effects of not only fine 
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structure which results in the sphtting of energies corresponding to various j-values, but also 
the hyper-fine structure and Lamb shift. The hyper- fine structure is due to the interaction 
between the spins of nucleus (proton in the hydrogen case) and electron, while the Lamb shift 
is caused by the virtual photon effects. To extract the necessary information about the fine 
structure, one can either carefully carry out a theoretical estimate to get rid of the influences 
from other effects or choose special transition modes where only the physical mechanism we 
concern plays a dominant role. In fact, we may ignore the effects of the hyperfine structure, 
because it is suppressed by a 



by the spin-orbit coupling 



'actor me/rrip ~ 1/2000 compared to the fine structure caused 



The energy gap between 2p-i and 2pi may be a good candidate to study the 



splitting 



18|. The 



id- 



induced by spin-orbit coupling, since they are in p wave and the Lamb shift is small 
level- crossing value between 2n3 and 2»i has been measured as 10969.20 ± 0.06MHz 

2 2 

Using this value we obtain a = -(8.93 ± 0.02) x 10"^ { a/b = -0.000893) which manifests 
the restoration degree of the SO (4) symmetry. 
There are three comments, 

(1) Since a is a very small number, Eq.(l6])is very close to Eq.(jl]), i.e. the symmetry 
breaking in the QED case is almost complete. 

(2) Especially, a is negative. Since this value is obtained by fitting data, we should suppose 
that it is valid. Question is what the negative sign implies? As we discuss above, the SO (4) 
symmetry demands a = 6=l/2, ifa were a small positive value, we could allege that the 
SO (4) symmetry were partly restored, but now it is a small negative value, what conclusion 
can we draw? We will come to this question at the last section. 

(3) Even though there is a scalar potential in Eq.([6]), the mechanisms which induce the 
Lamb shift and hyper-fine-structure are not involved in the Dirac equation (Eq.Q) as afore- 
mentioned. 



IV. ORIGIN OF SCALAR POTENTIAL IN QUANTUM FIELD THEORY 

In principle, we can make a general discussion about how the scalar potential emerges 
in the relativistic quantum mechanics. The least action principle in quantum mechanics 
tells us how to construct the Hamiltonian for a charged particle in electromagnetic field by 
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TABLE II: 



n 


n' 


3 




/ 


spectrum 


eigen-values (eV) 


1 





1 

2 


1 







13.605434057 


2 





3 
2 


2 


1 




3.400994431 


2 


1 


1 
2 


-1 


1 


2pi 


3.401039795 


2 


1 


1 
2 


1 





2s 1 

2 


3.401039795 



transforming operator p into p — eA and E into E — e(f){r) where 0(r) is the zeroth component 
of A'^^. The corresponding Dirac equation would be of the form [a ■ p + mj3 + (f){r)]tp = Eip 
where only the vector potential exists, because the energy E has the same Lorentz structure 
as (f) and it is exactly Vv{r) we concern. If the non-relativistic Hamilton can be written 



as H = ^J-p^ + (m + y(r))^, after the transformation, we can obtain a different equation 
[a ■ p + (m + V{y))l5\i\} = Etp, because now the potential V has the same Lorentz structure 
as the mass and then V is a. scalar potential. However there is no reason to write the 
Hamiltonian in that way. 

Now let us turn to the QED theory which is viewed as the correct theory to be responsible 
for all the electromagnetic processes, at the quantum field theory level. As aforementioned, 
the QED is a vector coupling theory, so that at the leading order the scalar potential does not 
exist. One can ask if the scalar potential emerges at higher orders. The answer is positive: 
actually we can deduce a scalar coupling from high order QED. Due to the fundamental 
coupling between electron and photon tp^'-fntpgA^, we achieve the vector potential from a 



single photon exchange 



20| 



When the high order effects are taken into account, namely, as we consider the loop 
contributions, the scalar coupling appears. The next to leading order Feynman diagrams 
with two photons exchanged between electrons are shown in Fig. 2 and 3. Omitting irrelevant 
factors in the amplitude, we can write it as 



El+J^ U(p)^U(P,h ^^±^7 f/rP)g^-g^- - 



^ Here to avoid unnecessary ambiguity, we do not use the terminology commonly used in electrodynamics: 
scalar and vector potentials which obviously have different meaning from what we refer in this work. 



8 



7a. 



FIG. 1: 

7m , Pi 7p _ 



Pb 



P2 lul 



The term 



FIG. 2: 



1 



can be rewritten as 

{Pi - m\){Pl - ml 
A simple relation 



'Pi -ml 



{Pa-PinPl 



-U{Pb)lpl^.U{Pa)U{Pd)i^^,U{P,)gP^g 



{Pa-PinPl 



1/ N 1/ N • 

7a 7/3 = 2 Wa7/3 " 7/37a) + 2 ( 7a 7/3 + 7/37a) = fi-a/J " ^C^a/?, 



In n ip 




Iv P2 luj 



Pd 



FIG. 3: 
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leads to 

mim2 



TTln 



-U{P,)gp^U{Pa)U{Pd)g^,U{P,)gf'^g 



-U{P,)U{Pa)U{P,)U{P^)- — — — , (12) 



i^Pf + ml){Pl + miy- "'^^ "^^^ '^^^^ ^^(P^-P,y^P,-P,)2-- 
which is exphcitly a 1 x 1 type couphng and eventually induces a scalar potential. The 
Feynman diagram in Fig. 3 also results in similar scalar coupling. 

Now let us turn to study the form of the induced potential. In fact, the integration 
is complicated and usually cannot be analytically derived. This is similar to deriving a 
potential between two nucleons induced by two-pion exchange. Instead of making a full 
integration over q which is supposed to correspond to the one-loop radiative correction to 
the scattering amplitude, we would like to derive the main piece of an effective potential 
between two charged particles. Let us focus on the two photon propagators in Eq. (fT2l) . we 
write the term as 

^ (13) 



g2(g + P)2' 

where q = Pi — P^ and P = Pf, — Pa, and the two fermion propagators are quenched 



The standard way 16|, ll7|] of obtaining a potential in QFT is to carry out a Fourier 



transformation on P in the transition amplitude. Now let us see the integration 



(2t)4 7 ,2|(,„ + P„)2-(q + : 

^ -A^V, (14) 



where the loop integration over four-momentum q can be carried out in the regular way. 
For the Fourier integration, following the standard, we set Pq = 0. Shifting the integration 
variable P into P' = P + q, we have 

dn'^ r pj(P'-q)r 



\d'P'. (15) 



W J [ql- 

The contour integration is standard and one immediately obtains the results in terms of the 
residue principle as 

(27r)4 g2 ~ 
10 



1 ( »-*-.,-„.,M. (16) 



r j (27r)4 g2 

It shows that the leading term in the induced potential is a generalized Coulomb-type, 
even though the numerator in Eq. (fT6ll is not a simple rational number, but an integration. 
It is worth noting that the numerator in Eq. ( fT6|) corresponds to the loop integration which is 
not Lorentz invariant because we have already chosen the rest frame of the the two-fermion 
system as our reference frame. To really calculate it, one may need to deal with possible 
infrared divergence which appears for soft photon propagation, but anyhow, it is a number. It 
is beyond the scope of this work, here we only show the characteristic of the potential induced 
by the two-photon exchange. Indeed, such a potential induced by the two-photon exchange 
is very complicated and it is in analog to the interaction between nucleons induced by two- 
pion exchange which has been thoroughly studied by many authors j2l|. The potential form 
is not easy to be framed, generally can be better described by the Bethe-Salpeter equation 



221. In Ref. 



22l |. the authors gave a very complicated potential form, but one can notice 
that the leading piece is still the Yukawa-type one which corresponds to the Coulomb-type 
in the case of two photon-exchange. 

As carrying out above theoretical calculations, we only consider the scalar Coulomb type 
potential in Eq. ffTTj) . but definitely when fitting data, contributions from all other pieces 
are involved altogether. It seems that our theoretical calculation is not complete. However, 
the contributions from other pieces are either of the same order of magnitude or, most 
probably, smaller than the contributions of the Coulomb piece. Therefore, as we discuss the 
symmetry breaking and restoration in this work, we only concern the order of magnitude, 
thus our results which only involve contribution from the Coulomb piece do make sense. 
In other words, our derivation is not rigorous, but is illustrative and persuasive. We will 
carry out a complete derivation of the two-photon induced potential in our later work. It 
also interesting to note that the sign of the Coulomb-type potential induced by two-photon 
exchange is opposite to that induced by one-photon exchange. It can naturally explain the 
minus sign of a obtained by fitting data. 

Because the scaler coupling arise at the next to leading order so the scalar potential 
should generally undergo a loop suppression. By a simple analysis of order of magnitude, this 
suppression factor is proportional to a/vr = ~ 0.00232 relative to the vector potential. 
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Definitely the concrete suppression factor sliould be derived by explicitly calculating the 
loop diagrams and numerically evaluate the contributions. But as the order O(a^) and even 
higher orders can contribute to the scalar potential, we would like to obtain the suppression 
factor a/h phenomenologically. The number on a/6 which is obtained by fitting the spectrum 
of hydrogen and the value was given in last section, one can notice that it is consistent with 
our qualitative analysis. 



V. CONCLUSION AND DISCUSSION 



It is a well-known fact that symmetry plays an important role in physics. Not only the 
geometric symmetries are obviously important, but also dynamical symmetries are funda- 
mental and crucial for determining the physics. For example. It has been a long time, people 



notice that there is an extra symmetry [R, H] =Q with R = 2^(pxL — Lxp) — for the 
hydrogen atom. This symmetry finally leads to the degeneracy of the hydrogen levels 23|] . 
The symmetry corresponds to two "angular momenta" and the generators are respectively 

2 ' 

and 

2 ' 

where K = ^/^^R. The two "angular momenta" symmetries are involved in a larger group 
S0(4). 

Similarly, it is indicated that the Hamiltonian of harmonic oscillator has an SU(3) dy- 
namical symmetry. Since the Schrodinger equation is the relativistic approximation of the 
relativistic Dirac equation for fermions, one can ask if the larger symmetries are maintained 
in the corresponding Dirac equations? The symmetry causes degeneracy and the symmetry 
breaking leads to splitting of energy levels, thus to answer the question, we need to study 
the observational energy level splitting. Ginocchio 6, 7| and Zhang et al. indicate that to 
maintain the higher symmetry in the Dirac equations for harmonic oscillator and hydrogen 
atom, the scalar and vector potentials must have the same weight, i.e. the combination is of 
the form (1 + (3)V{r). If one changes the combination into (a + h(3)V{r) with a + 6 = 1, as 
he reduces the Dirac equation into the non-relativistic Schrodinger equation, at the leading 



12 



, / , 

order, the combination does not induce distinction because /? = , only the left 

-/ 

corner entry plays a role. But to the next order, the right corner entry joins the game and 
a and h terms would have different contributions and the superficial symmetry is broken 
unless a = h = 1/2. If so, the spin-orbit coupling does not exist and fine structure (energy 
splitting) disappears. It obviously contradicts to the experimental observation. 

The underlying theory is QED which is a vector-like theory, and its leading order contri- 
bution to the potential is via the one-photon exchange which definitely results in a vector 
potential in the Dirac equation. Even though the higher order QED effects would induce the 
scalar potential, it should be much suppressed comparing to the leading term. The data tell 
us that the fraction of the scalar potential is at order of 0.1% which is consistent with our 
estimation based on the quantum field theory. We obtain the value of a as a = — (8.93±0.02) 
( a/h = —0.00089) by fitting the hydrogen spectrum, and it seems that the value further 
deviates from a = b = 1/2 which is demanded by the S0(4) symmetry. Does the negative 
sign really mean that the symmetry is not partly restored? Since the value of a is too small, 
some theoretical and experimental uncertainties can cause a small fluctuation in the value, 
so that we cannot yet draw a definite conclusion that the symmetry is not restored. To 
confirm its sign, one must pursue this study in a case where the restoration might be sizable, 
and no doubt it is the QCD case. 

Since the higher order contributions of QCD would be much more important because a.^ 



is large at lower energies, so one can expect that in t 



might be significant. Castro and Franklin 



24, 



25 



le QCD case, such symmetry restoration 



26 



271] have worked out the Dirac equation 



where the authors considered an arbitrary combination of the scalar and vector potentials 
for both the Coulomb and linear pieces and obtain the wavefunction of ground state of a 
qq bound state. Interesting aspect is that they find that the dominant part in the linear 
confinement potential is the scalar part. That conclusion is obviously different from the 
QED case. Since the linear potential is fully due to the long-distance effects of QCD and 
may be caused by the non-perturbative contributions, we cannot derive it from QFT so far 
(in the the lattice theory, the area theorem may imply a linear potential form, but cannot 
distinguish the vector and scalar fractions at all.), the argument about the origin of scalar 
potential given in last section becomes incomplete. We will further pursue this interesting 
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issue and it is the goal of our next work 28 ] . 



Our conclusion is whether the high symmetries in the Schrodiner equation may exist in 
the corresponding Dirac equation needs careful studies. In the case of hydrogen atom, it 
definitely does not. We show that the symmetry breaking is almost complete, from the 
experimental and theoretical aspects. 

Indeed it would an interesting problem because we explicitly show the breaking for hydro- 
gen atom, but for other kinds of potentials which may be induced by various mechanisms, 
especially the phenomenological potentials such as the Cornell-type potential which is respon- 
sible for the quark confinement, whether such symmetry can be kept in the Dirac equation 
is worth further investigation. 
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